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ABSTRACT 


The  Spectral  Deconvolution  Analysis  Tool  (SDAT)  software  was  developed  to  improve  counting  statistics 
and  detection  limits  for  nuclear  explosion  radionuclide  measurements.  SDAT  utilizes  spectral 
deconvolution  spectroscopy  techniques  and  can  analyze  both  p-y  coincidence  spectra  for  radioxenon 
isotopes  and  high-resolution  HPGe  spectra  from  aerosol  monitors. 

The  SDAT  tool  has  been  integrated  into  a  standalone  graphical  user  interface  (GUI).  This  GUI  may  import 
spectra,  analyze  the  data  for  Xe  concentrations,  and  graphically  display  the  results.  This  tool  has  been 
tested  with  data  generated  via  MCNPX  models  as  well  as  real  data. 
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OBJECTIVES 

Environmental  xenon  sampling  and  measurement  units  are  used  in  nuclear  weapons  test  monitoring 
networks  because  radioxenon  isotopes  may  be  all  that  is  released  from  tests  performed  deep  underground  or 
underwater.  Some  environmental  xenon  sampling  and  measurement  units,  like  the  Automated  Radioxenon 
Sampler-Analyzer  (ARSA)  developed  by  Pacific  Northwest  National  Laboratory  (PNNL)  and  the  Swedish 
Automated  Noble  Gas  Unit  (SAUNA)  developed  by  the  Swedish  National  Defense  Research  Establishment 
(FOI),  use  p- /coincidence  detectors  that  are  energy  dispersive  on  both  the  /?and  /energy  axes.  Signals 
from  four  radioxenon  isotopes  ('  lmXe,  133Xe,  1  mXe,  135Xe)  comprise  a  sample  spectrum.  Under  poor 
operating  conditions,  a  few  radon  daughters  (214Pb,  2l4Bi)  may  interfere  with  the  sample  spectrum. 

Applying  conventional  region-of-interest  (ROI)  spectrum  analysis  algorithms  to  such  3-D  spectra  results  in 
relatively  high  minimum  detectable  concentrations  (MDCs)  due  to  the  subtractive  process  of  determining 
net  counts  in  the  ROI. 

To  overcome  this  problem,  a  project  team  from  The  University  of  Texas  in  Austin  is  developing  software 
to  deconvolve  the  3-D  sample  spectra  into  the  most  probable  combination  of  signals  using  the  non-negative 
least-squares  method.  This  method  can  use  the  entire  signal  from  each  radionuclide  and  consequently  may 
improve  the  signal  to  noise  ratio.  It  has  been  postulated  that  the  use  of  such  an  algorithm  will  result  in  a 
better  ability  to  resolve  spectral  interferences  and  thus  improve  counting  statistics  and  detection  limits  for 
nuclear  explosion  radionuclide  measurements.  The  SDAT  GUI,  expands  upon  a  previously  developed 
algorithm,  the  Multiple  Isotope  Component  Analysis  (MICA)  tool. 

The  SDAT  tool  was  tested  against  the  CO  RIANT  code  (Foltz  Biegalski,  2001)  on  data  produced  from  an 
automated  radioxenon  sampler/analyzer  (ARSA)  unit  when  it  was  stationed  in  Guangzhou,  China.  The 
CORIANT  tool  utilizes  a  region  of  interest  approach  to  the  spectral  analysis  while  SDAT  utilizes  a  spectral 
deconvolution  approach.  Samples  were  selected  that  contain  clear  radioxenon  signals.  While  a  true  test  of 
SDAT  will  be  the  analysis  of  samples  with  low  count  rates,  this  test  was  conducted  as  a  first  step  to  show 
that  the  two  methods  produce  similar  results. 


RESEARCH  ACCOMPLISHED 

This  occurred  in  three  steps.  The  first  step  was  the  creation  of  the  SDAT  tool.  The  second  step  was  the 
creation  of  standard  library  inputs  for  the  SDAT  code,  which  was  accomplished  via  MCNP  simulations. 
The  third  step  was  to  analyze  samples  with  SDAT  and  compare  the  results  to  those  produced  with 
CORIANT. 

Step  1:  Creation  of  SDAT  GUI 

The  detector  response  for  a  sample  consists  of  a  255x255  matrix  of  numbers.  Each  entry  represents  the 
number  of  counts  registered  in  a  certain  E/  Ep  bin.  Each  row  represents  a  /-channel  bin  and  each  column 
represents  a  /^-channel  bin.  We  will  refer  to  this  structure  as  a  histogram.  Other  information  contained  in 
the  sample  file  includes  calibration  information  and  other  sample  characteristics,  e.g.,  the  sample  Xe  gas 
volume  from  which  the  total  sampled  atmospheric  volume  is  calculated. 

The  sample  histogram  will  be  deconvolved  into  individual  isotopic  responses  using  the  SDAT  concept  to 
determine  atmospheric  activity  concentrations  for  each  radioxenon  of  interest.  To  do  this,  however,  we 
must  have  calibrated  histograms  of  all  the  possible  individual  signals  that  can  make  up  a  sample  histogram. 
These  histograms  should  have  the  same  size  and  calibration  characteristics  as  the  sample  histogram  in 
addition  to  good  counting  statistics.  Therefore,  we  need  the  following  detector  response  matrices  with  their 
associated  activities: 

1  33m 


133  y  135  y 

^1255x255,  1  ^1255x255 ,  and  ^ 


1255x255 


131m  y 

^^255x255  >  - 
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Each  detector-response  histogram  can  be  generated  using  a  detector  modeling  program  like  MCNP  or 
acquired  by  counting  a  calibration  source  on  the  detector.  Since  the  ARSA  has  four  beta  detector  cells,  the 
above  detector  response  histograms  would  need  to  be  generated  or  acquired  for  each  beta  cell.  In  addition, 
the  energy,  resolution,  and  efficiency  calibrations  have  to  be  determined  for  each  beta  cell.  With  all  of  this 
information  at  hand,  and  knowing  in  which  detector  cell  the  sample  was  counted,  the  following  algorithm 
can  be  applied  to  a  sample  for  determining  the  atmospheric  radioxenon  concentrations. 

The  SDAT  code  was  developed  to  input  the  library  and  sample  files.  These  files  are  text  files.  The  current 
version  of  SDAT  just  reads  the  sample  spectrum  and  does  not  read  any  header  information  that  may  be 
associated  with  Sample  Pulse  Height  Data  (SPHD)  formatted  files.  The  open  sample  window  is  shown  in 
Figure  1. 

Once  the  proper  files  have  been  selected,  SDAT  analyzes  the  data  and  displays  the  results.  Figure  2  shows 
the  results  display.  The  sample  is  displayed  as  a  2-D  plot  in  the  top  left  corner.  Below  the  sample  plot  are 
three  tables.  The  top  left  table  is  the  Axis  Parameters  table  which  allows  users  to  modify  the  minimum  and 
maximum  x  and  y  intercepts  of  all  the  plots  in  the  SDAT  window.  By  default,  they  are  set  at  0  and  250.  The 
Parameter/Value  Table  is  situated  to  the  right  of  the  Axis  Parameters  table  and  contains  a  summary  of  the 
parameters  specified  in  the  Open  Single  window,  i.e.,  the  sensitivity  number,  the  sample  file  name,  and  the 
library  file  names  for  each  radioxenon  of  interest.  This  table  is  located  in  the  bottom  left  corner  of  the 
SDAT  window.  It  contains  the  concentration  coefficients  and  errors  of  each  radioxenon  of  interest 
calculated  with  and  without  the  use  of  data  weighting  via  the  weighting  matrix  described  earlier.  The 
resulting  multiplier  coefficients  can  be  used  to  determine  the  activity  of  each  radioxenon  using  equation  1 . 


c, 


M,A> 

V 


(1) 


where 

C,  =  Sample  activity  concentration  of  radioxenon  isotope  /; 

Mj  =  Multiplier  coefficient  of  radioxenon  isotope  i  as  determined  by  SDAT ; 
A  j  =  Activity  represented  by  library  file  of  radioxenon  isotope  /;  and 
V  =  Air  volume  sampled. 


The  activity  concentration  uncertainties  of  each  radioxenon  can  be  determined  using  equation  2. 


a 


c, 


(2) 


where  <JV  is  the  uncertainty  of  variable  x.  All  other  variables  have  been  defined  above. 


The  weighted  result  and  unweighted  result  plots  are  located  at  the  top  middle  and  top  right  corner  of  the 
SDAT  window.  The  plots  are  calculated  as  shown  in  equation  3. 

(!£ [lib], ,)[#],  ,3, 

where  C,  are  the  concentration  coefficients  for  radioxenon  isotope  i,  [ lib ],  are  the  library  files  for 
radioxenon  isotope  i,  and  [cj  is  the  weighting  matrix  for  the  weighted  result  plot  and  the  identity  matrix  for 
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the  unweighted  result  plot.  These  plots  may  be  compared  with  one  another  and  the  original  sample  plot  to 
see  how  well  the  SDAT  calculated  results  model  the  sample  data. 

The  Weighted  Residual  and  Unweighted  Residual  Plots  are  located  directly  below  the  Weighted  Result  and 
Unweighted  Result  Plots.  The  plots  are  calculated  shown  in  equation  4. 

[Sample]- ((^.C,  [ lib\  )•  [£])  | ,  (4) 

where  [Sample]  is  the  original  sample  data  and  all  other  symbols  have  been  previously  defined.  These  plots 
show  the  differences  between  the  original  sample  data  and  the  calculated  models  of  the  sample  as  shown  in 
the  weighted  result  and  unweighted  result  plots.  SDAT  also  allows  for  sample  re-analysis  may  with 
different  input  data.  Results  may  also  be  saved. 

Over  time,  the  energy  vs.  channel  calibration  of  both  beta  and  gamma  detectors  can  change.  If  the  changes 
are  significant,  the  SDAT  program  will  return  large  errors  because  the  sample  signal  will  no  longer  match 
those  in  the  library  files.  In  such  cases,  the  sample  spectra  should  be  rescaled  to  match  the  energy 
calibration  properties  of  the  original  library  files.  Alternatively,  the  library  files  can  be  rescaled  to  match 
the  sample  calibration  properties  and  saved.  This  program  can  also  be  used  to  compress  or  expand  the  data 
as  desired.  Figure  3  shows  the  SDAT  rescaling  window. 

Step  2:  Creation  of  Xe  Library  Files  via  MCNP 

Library  files  were  created  for  l31mXe,  133mXe,  133Xe,  and  135Xe  via  a  combination  of  MCNPX  and 
MATLAB  post  processing.  This  is  the  same  method  as  reported  in  Haas  et  al.  (2007).  The  library  files  were 
created  for  each  of  the  four  cells  in  the  (3-y  coincidence  detector  system  utilized  in  the  ARSA  system.  The 
library  files  were  then  adjusted  to  match  the  energy  calibration  of  the  samples  collected  in  the  field.  Figure 
4  shows  a  comparison  between  the  synthetically  generated  133Xe  library  file  and  a  sample  with  L33Xe 
collected  by  the  ARSA  in  Guangzhou,  China. 


Step  3:  Comparison  of  SDAT  with  CORIANT  Results 

As  an  initial  quality  control  check,  five  spectra  were  analyzed  wit  SDAT.  These  samples  collected  with  an 
ARSA  unit  in  Guangzhou,  China  were  initially  analyzed  by  CORIANT  at  the  Center  for  Monitoring 
Research  in  Arlington,  VA.  Figure  5  illustrates  the  comparison  that  was  made  between  the  SDAT  results 
and  the  CORIANT  results.  The  133Xe  results  compare  very  well.  Since  these  results  are  from  samples 
collected  in  the  field,  the  correct  quantity  of  133Xe  in  the  sample  is  not  known.  This  just  serves  as  a  check  to 
make  sure  that  the  SDAT  tool  produces  acceptable  results  for  samples  with  strong  Xe  signals. 

The  analysis  of  131mXe,  133mXe,  and  135Xe  were  also  performed.  These  results  did  not  compare  as  well  as  the 
133Xe  results  shown  in  figure  5.  This  issues  with  the  131mXe  and  133mXe  was  primarily  a  detection  issue. 
CORIANT  reported  131mXe  and  133mXe  present  in  all  of  the  samples  analyzed.  However,  SDAT  calculated 
131mXe  to  be  present  in  four  out  of  the  five  samples  and  did  not  calculate  133mXe  to  be  present  in  any  of  the 
samples.  Discrepancies  in  these  results  are  largely  attributed  to  having  correct  (3  channel  calibrations,  which 
are  questionable  for  both  the  CORIANT  and  the  SDAT  analysis.  Visual  inspection  of  the  data  does  not 
show  any  clear  131mXe  or  133mXe  signals. 

Another  issue  raised  by  this  comparison  is  the  error  calculation  performed  in  SDAT.  The  error  bars 
produced  for  this  comparison  appear  to  be  smaller  than  expected.  While  the  error  calculations  may  be 
correct,  efforts  should  be  made  to  review  the  error  calculation  in  SDAT. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  first  version  of  the  SDAT  code  was  completed.  Samples  were  analyzed  with  the  code  using 
MCNPX/MATLAB  generated  library  files.  Results  from  this  comparison  agreed  with  CORIANT  results 
from  the  Center  for  Monitoring  Research.  Future  work  for  the  SDAT  code  will  include  an  incorporation  of 
activity  calculations  as  well  as  an  ability  to  read  SPHD  file  headers. 

The  highest  priority  work  will  entail  the  creation  of  radioxenon  isotopes  via  The  University  of  Texas  at 
Austin  TRIGA  reactor.  These  spectra  will  be  important  to  use  as  library  files  for  SDAT.  They  will  also 
allow  for  advanced  testing  of  SDAT  with  low  counting  statistics  samples.  Since  the  xenon  concentrations 
may  be  experimentally  controlled,  the  accuracy  and  precision  of  SDAT  may  be  determined. 
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Figure  1.  SDAT  GUI  open  single  window. 
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Figure  2.  SDAT  GUI  sample  analysis  results  display. 
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Figure  3.  Display  after  rescaling  a  sample. 
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Channel 

Figure  4.  Comparison  of  |3  and  y  compressed  spectra  from  MCNPX/MATLAB  generated  133Xe 
library  file  and  from  sample  collected  in  Guangzhou,  China. 
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Figure  5.  Comparison  of  results  between  SDAT  and  CORIANT  for  five  samples  from  the  ARSA  in 
Guangzhou,  China. 
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